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We study the angular momentum mixing effects in the color superconductor with non-spherical 
pairing. We first clarify the concept of the angular momentum mixing with a toy model for non- 
relativistic and spinless fermions. Then we derive the gap equation for the polar phase of dense 
QCD by minimizing the CJT free energy. The solution of the gap equation consists of all angular 
momentum partial waves of odd parity. The corresponding free energy is found to be lower than 
that reported in the literature with p-wave only. 
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^ , I. INTRODUCTION 

. The properties of quark matter at extreme conditions have been an active research area both theoretically and 
^ ' experimentally. At high temperature, the quark-gluon plasma(QGP) has long been searched by colliding two nuclei 
at sufficiently high energy. On the other side, we expect that quark matter becomes color superconducting through 
a phase transition at high baryon density but low temperature[l|, [3, 0, 0j S [1] i which is the typical condition inside 
compact stars. 

^ • In a typical metallic superconductor, the electrons pair with equal chemical potential near the Fermi surface. 
1^ ' The situation with a quark matter, however, is much more complicated. While the quark-quark interaction favors 
> pairing between quarks of different flavors, the mass difference among u, d and s together with the charge neutrality 
requirement induces a substantial mismatch among their Fermi momenta at the baryon density inside a compact star 
, and thereby reduces the available phase space for Cooper pairing. A number of exotic color superconductivity phases 
^ in the presence of mismatch have been proposed in the literature B, IC , Hi. b ut a consensus point of view of the 
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true ground state has not been reached. The single flavor pairing which is free from the Fermi 

0^ momentum mismatch, is an interesting alternative in this circumstance and will be considered here. Since the quark- 



quark interaction is attractive in the color anti-triplet channel, the color wave function of the pair is anti-symmetric. 
For the equal helicity pairing to be considered in this article, the parity of the orbital wave function has to be odd 
as required by the Pauli principle. Except for the color-spin-lock phase examined in jl^, the energy gap will not be 
spherical. The odd parity prevents the diquark wave function from realizing the full pairing potential. The energy 
scale of the color superconductivity is therefore reduced. 

At ultra-high baryon densities, asymptotic freedom of QCD ensures the validity of the weak coupling expansion, 
which has been carried out for CSC by a number of authors^ [13, [H, [MUS [HI, Hi [H, [H, [H, [ll] . The dominant 
pairing interaction is mediated by one-gluon exchange and can be decomposed into partial waves, as is shown in 
Eq. (|59p below. Quantitative results of the transition temperature in the equal helicity channel with an arbitrary 



Ti^) = 5I2.^( exp [ - ^ + 7 - ^(vr^ + 4) + 3cj\ (1) 

where Nf is the number of flavor, N is the number of colors, fj, is the chemical potential, g is the running coupling 
constant of QCD and 7(= 0.5772...) is the Euler constant. The J-dependent constant 

0, for J = 0, 

CJ={ J (2) 
-2^1 for J > 0. 

n=l 
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The equal helicity pairing of odd parity picks up Tc^^ as the transition temperature. We have Tc^^ = e ^TjP^ ~ 
2.479 X 10-3ri°^ Natural analo gy is drawn with the superfluidity of He'^. But important difference between the 
pairing potential in quark matter and that in He'^ has to be considered before ascertaining the angular dependence of 
the energy gap. The forward singularity of the one-gluon exchange renders the pairing strength equal for all partial 
waves (the same leading order term inside the bracket of ([T])). The paring potential in He^, however, is entirely in the 
channel of J = 1. 

The transition temperature was determined from the pairing instability of the diquark scattering amplitude in the 
normal phase. In a perturbative treatment, the Dyson-Schwinger equation for the scattering amplitude is linear and 
the partial wave decomposition in Eq.(IT]) is legitimate. This is not the case with the gap equation below Tc, which 
is nonlinear. A non-spherical solution of the gap equation, in general, will be a mixture of different partial waves 
unless the the pairing is dominated in one angular momentum channel. The adjective "non-spherical" refers to the 
magnitude of the energy gap of the quasi-particle. So the s-wave and the CSL gaps are spherical and are therefore 
free from the angular momentum mixing. The gap equations in these cases are linear with respect to the angular 
dependence. A nonsphcrical gap function, (pM , for the equal helicity pairing at T < Tc contains all spherical harmonics 
of odd parity with the same azimuthal quantum number M . We have 

(^M = 02Sc/m(p) (3) 

where 02SC is the gap function of 2SC in the absence of the mismatch, p the direction of the relative momentum of 
the two quarks in a Cooper pair and the angular dependent factor 

Im{p)^ E bjYjMip) (4) 

J=l,3,5,... 

with J the total angular momentum of the Cooper pair. Carrying the formulation of He"^ over to QCD amounts to 
drop all higher multipoles except that of J = 1, which will not satisfy the gap equation of QCD. It was argued in 
the literature that hi — 0{1) but hj = 0{g) for J 7^ 1, This, as will be shown below, is not the case. Instead, we 
find that the function /m(p) satisfies a nontrivial integral equation and thus bj = 0(1) for all odd J's. Therefore the 
angular momentum mixing does occurs in the subleading order of the gap function. The angular momentum mixing 
will modify all non-spherical "spin-1" CSC examined in the literature, we shall focus our attention in this paper to the 
equal helicity pairing with zero azimuthal quantum number, i.e. the analog of the polar phase of He''. The subscript 
M of (pM and /m will be suppressed below. Even though this phase is unstable, it is the simplest one to illustrate 
the mixing mechanism. 

The current work is organized as follows. In the next section, we shall clarify the concept of the angular momentum 
mixing with a toy model of non-relativistic and spinless fermions. In the Sect. Ill the gap equation for the single flavor 
CSC will be derived by minimizing the CJT free energy of QCD. This gap equation will be reduced to an nonlinear 
integral equation for the function /(p) in the subsequent section and the numerical solution will be presented in the 
Sect, V. We conclude the paper in the Sect. VI. Some technical details are deferred to the Appendices. Our units 
are h = c = ks = ^ and 4- vectors are denoted by capital letters, K = = {ko,k) with ko the Matsubara energy, 
which becomes continuous at Tc = 0. Throughout the article, we shall follow the deflnition of the leading order and 
the subleading order in [6] . Upon taking the logarithm of the transition temperature or the magnitude of the gap 
function, the order 0{-) will be referred to as the leading one and the 0(1) term to the subleading one. 



II. A TOY MODEL WITH ANGULAR MOMENTUM MIXING 

To clarify the concept of the angular momentum mixing, we consider a toy model of nonrelativistic and spinless 
fermions. In terms of the creation and annihilation operators, the model Hamiltonian reads 

H = Y^ epuiaj; - -^Yl ■ p')4"-p"-P'"p' 
p p.p' 

2 

where tp = ^ — /i with m the mass and /i the chemical potential, A > is a coupling constant, fl is the normalization 
volume and the summation '}2'pp' extends to states with |ep| < ujd and |ep/| < ujd with tj/j a UV cutoff (Debye 
frequency for electronic superconductors). The angular dependent form factor V{p ■ p') can be expanded in series of 
the Legendre polynomials, 

V{p-p')^J2'^2J+l)vjPj{p-p'). (6) 

,7=0 
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Introducing the order parameter 

X{p) =< la-pUffl > (7) 

with I > the ground state and expanding the interaction term of ([5]) to the Uncar order of the fluctuation a-.pap— < 
\a^pap\ >, we end up with the mean-field Hamiltonian 

p,\ep\<uJo P p,\ep\<ujn 

where we have introduced the gap function via 

<^(p^ = ^ E yip-i^)xm- (9) 

p'.\ep,\<uJo 

We have xi^P) — ^xip) ^^'^ 4'{~P) — ~4>{p) foUowing from their definitions. Upon a Bogohubov transformation, we 
find that 

m - (10) 



with Sp = Jtp + and the ground state energy 



So = ep + ilF (11) 

p,ep<0 



with F the condensation energy density given by 



F 



2 .. z,c 

p,\ep\<ujD ^ 



Substituting pO|l into ([9]), we obtain the gap equation 



p',\tpl\<u]D 



3 1 

In the weak couphng hmit, ujjj « fi and XDp << 1 with Dp ~ "^^^ the density of states on the Fermi surface, 
but the magnitude of V{p ■ p') remains of order one. We have 



p,\<ip\<ulD 



Also, the support of the gap function extends only to a narrow band around the Fermi surface. We may ignore the 
dependence of </) on the magnitude p = |p| and switch the argument of (f) from p to p. Following (|14p , the integrations 
over p in (|12p and (jl3[) can be carried out readily and we end up with 

and 



^ip) = ^ y '^ p'^^p ■ p>^p'^ ^""WFYl' ^^^^ 

The gap equation pH]) is nonlinear because of the logarithm. 

In what follows, we consider two extremes of V{p ■ p'), each of which gives rise to an exact solution to the gap 
equation (jl6p . We present only the solution that is invariant under time reversal, i.e. the one with zero azimuthal 
quantum number. 
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Case I: 

V{p-p')^W,ip-p')^3p-p'. (17) 

It corresponds to the partial wave expansion ([6]) with vj = 1 and vj — for J ^ 1. The angular dependence of the 
pairing force in He'^ is of this type. The gap equation reads 

and its solution of zero azimuthal quantum number is given by 

(j){p) = 00 cos e = (j)oPi (cos 9) (19) 
with 9 the angle with respect to a prefixed direction in space and 

(j)o ^2ujDe~^~^K (20) 

The condensation energy density 

F = ^T^e ^"i-^^ ~ -0.6492w|,i:ii.e . (21) 

o 

This solution corresponds to the polar phase of He'^ "23?!. Since the gap function contains only the partial wave of 
J = 1, there is no angular momentum mixing. The additional term in the exponent of (|20p . i, comes from the 
logarithm of (fT8)) . 
Case II: 

V{p ■ p') = 4n5\p - p') = ^(2 J + l)Pj{p ■ p'). (22) 

J 

This corresponds to a singularity of the two body scattering amplitudes in the forward direction. We have vj — \ for 
all J in ([ni . The last step of follows from the addition theorem and the completeness of the spherical harmonics. 
The gap equation (fTH]) becomes 

\Df 2ujd 

which imphes a constant and yields a solution of odd parity and zero azimuthal quantum number. 

^ ^ cxD ('2n ll!! 

m - 2a;z5e-^sign(cos^?) = 2u;De-^ ^(-l)"(4n + 3)-^^^— ^P2„+i(cos 0). (24) 

The condensation energy density in this case reads 

F=-ujlDpe^^. (25) 

We refer to this case as the case with the angular momentum mixing because the gap function ([24)) contains all partial 
waves. Carrying the solution of the case I to the case II amounts to drop all partial waves other than that of J = 1 
and would lead to a lower magnitude of the condensation energy (j2ip . 

The case with QCD is similar to the case II above since the forward singularity of the diquark scattering renders the 
pairing strength of all partial waves equal to the leading order. The running coupling constant g of QCD corresponds 
to A here and the angular momentum mixing shows up in the 0(1) term of ln|(/)|. Therefore we expect angular 
momentum mixing to the subleading order of the angular dependence of the gap function. Besides being an ultra 
relativistic system, the CSC of QCD differs from the toy model considered above in two aspects. The forward 
singularity of QCD also brings about the energy dependence of the gap, so the gap equation will be replaced by 
the Eliashberg equation derived by minimizing the CJT effective action of QCD. Secondly, the pairing strength of 
each partial wave does fall off with an increasing J in the sub-leading order of the pairing potential. It is this falling 
off that makes the amount of the angular momentum mixing numerically small for the solution considered in this 
article. 
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FIG. 1: The sunset diagram of Ea. (|29|l . Straight line and wavy hne denotes quark and gluon propagators respectively. 



III. DERIVATION OF THE GAP EQUATION FROM THE CJT FREE ENERGY 

The QCD Lagrangian for one flavor of massless quark is given by 

L = ^{ij'^Di^i + fi'jo)ip — ^G^'^G"^ + renormalization counterterms (26) 

where, tp is the quark spinor in Dirac and color space and -0 = ■0^70- The covariant derivative acting on the fermion 
field is Df^ = dfj_+igTaA^, where g is the running coupling constant, A° is the gauge potential, Ta = Xa/'2{a = 1, 8) 
is the a-th SU{3)c generator with Aq the a-th Gell-Mann matrix. = df,A'^ - d^A'^ + gf°-'"'A''^Al is the field 

strength tensor. Introducing the Nambu-Gorkov spinor 

*=fiV * = (Vi,^c) (27) 



where ipc = CiP^ is the charge-conjugate spinor and G = i'j^j^, the CJT effective action reads [gI Isoj 

r[D, S] = -i{TrlnL>-i + TiiD^^D - 1) - TrlnS^^ - Tr{So^S ~ 1) ~ 2r2[D, S]} (28) 

where D and S are the full gluon and quark propagators, Dq^ and Sq^ are the inverse tree-level propagators for 
gluons and quarks, respectively. r2 is the sum of all two-particle irreducible(2PI)vacuum diagrams built with D, S 
and the tree-level quark-gluon vertex T. We have 

T2^-^TT{DtSrS) + ..., (29) 

where the first term corresponds to the sunset diagram of Figlljand the contribution from ... is beyond the subleading 
order of the gap function[6] . 

The stationary points of the CJT effective action are determined by 

which gives rise to the Dyson-Schwinger equation for gluons and quarks, 

where 11 and E are the gluon and quark self-energy defined via = Do ^ + n and = S^^ + E. Instead of 
solving the two equations of ((3T|) simultaneously, we shall reduce the CJT effective free energy with the aid of the 
first equation, leaving the gap function arbitrary. The gap equation ( which is the Nambu-Gorkov off diagonal part 
of the second equation of (l3T|) ) will be derived after the CJT free energy is fully simplified under the weak coupling 
approximation. 

Substituting the first equation of ((3T|) into ((28| . the second term in Eq. ((28)) cancels the last term. We have 

r[D, S] = -^[TrlnD-^ - TrlnS'^^ - Tr(S'-^S' - 1)] (32) 
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In Nambu-Gorkov space, the inverse free quark propagator is 



[G+]-i 



where 



On writing the quark self-energy 



the full quark propagator, 



^0 -[ [G^]-i ' ^^^^ 



[Go±]-^ = (po±Ai)7o-7-P- (34) 



can be obtained explicitly by inverting the matrix ^ + E. 

For the single flavor pairing, the simplest choice of the off-diagonal block of Eq. ([55)1 reads 

$+(P) = i075A2 (37) 

and = $~(see Theorem 2 in[26j), where A2 is the 2nd Gell-Mann matrix and is a function of the energy and 
the momentum, i.e. (j) = 4>iPojP)- 4' is even in po and odd in p. By using the energy projectors of massless fermions 
Ap = (1 ± 7o7 • p)/2 and ignoring the contribution from the wave-function renormalization, the NG blocks of the 
propagator ([55]) take the form 

^ -Pl-{PT M)^ - <t>^\l ^" + - (p ± M)^ - 0^Ai 



" p2_(p±^)2_^2^2 P ' p2_(p^^)2_^2^2 P 

Because of the A2 of ([57]) . the excitation in the third color direction is ungapped. 

Now, we proceed to simplify the CJT free energy under the weak coupling approximation. Denote by r„ the free 
energy density of the normal phase, we have 

r = r„ + VLF (40) 

where the condensate energy density 

F = --L[Trlni:>-i - TrlnD-^ - TrlnS'-^ + TrlnSo"^ - 1:t{Sq^S - 1)] (41) 

is the part of V responsible to the gap equation. Following the procedure of (sil. [s^. we approximate 

TrlnZ?"^ - Trln[D^^] ~ Tt:[D,,5II\ (42) 

where 

m = n-n„, (43) 

with n„ the hard-dense-loop (HDL) resummed gluon self-energy in normal phase and Dn the corresponding HDL 
gluon propagator. In the Coulomb gauge, the HDL gluon propagator is 

D^^^o{K)^Di{K), Dr^MiK) ^ Dr^M = 0, D^^,, ^ [5,, - kk)Dt{K) (44) 

where Di^t are the longitudinal and transverse propagators respectively and are diagonal in adjoint color space, 
i.e. Df\ = (5°^_D;.f. Consequently, we only need the 00-component, H°"(iir), and the transverse projection of the 
ij-components, 

((5,y - hkjW^{K) = W{K) - hkjW^{K) (45) 
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The gluon self-energy in super phase reads 



p.P' 



where K = P — P' and 



r"" 
r'' 



(46) 



(47) 



with = 7^ra and = -l^Tj . Substituting Eq. ^ into Eq.dlSl), we find that Nambu-Gorkov space, 
KliK) = iT^{Tr[r;^G+(P)r;^G+(P')]+Tr[ff^G-(P)f|;G-(P')] 



p.P' 



-Tr[r^^S-(P)f,^S+(P')] + Tr[f^^S+(P)r,^S-(P')]} 



(48) 



Tr[rj^G+(P)r:;G+(P')] = Tr[7^A+7o7'^A+7o]u;+(P,P'), 



p.P' 



Since the HDL gluon propagators are diagonal in color space, we only need the diagonal terms of Eq. (|48|) to deal with 
Eq. (|42ll . The explicit form of each diagonal term of (|48p reads 



(49a) 
(49b) 
(49c) 
(49d) 



Tr[f^^G-(P)f:^G-(P')] = ^Tj] Tr[7^Ap-7o7^A;,7o]«^-(P,P'), 



p.P' 



Tr[r^^S-(P)f:^S+(P')] = -^T^ Tr[7^5A+7''75A;,]z«.(P,P'), 



p.P' 



T^[f^S+(P)r::S-(P')] = -^T 5] Tr[7^5A,-7^75A+>. (P,P'). 



p.P' 



where the repeated color indexes on LHS are not to be summed. The quantities and w on RHS of Eqs. (l49aj|49d|) 
are given by 



' 1 Po ± Cp Po ± ep/ 



< 4 



Po ± ep Po ± Cp' , Po ± ep p'q ± Cp' 



pI - ^1 p'o - 4' Po - 4 p'o - 4' 



,a = 1,2,3 



,a = 4,- • -,7 



(50a) 



and 



r 1 0(P)0(P') 



a = 8 



Wn = ( , a = 4, • • •, 7 

1 0(P)0(P') 



(50b) 



I 6(pg-e2)(p.2_^2,) 



,a = 8 



where ep = p — n and ep = \/ {p — [i)^ + (IP'{P). Since the dominant contributions in the weak coupling arise from the 
quasiparticles, we have ignored the contributions from the quasi- antiparticles in the calculations above. The trace 
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over Dirac space is straightforward 

Tr[70A±7o7*'A±7o] = -Tr[7075A±7°75AJ,] = 1 + p • p', (51a) 

J2 Tr[7*A±7of A±7o] = ^ Tr[7^75A±7'75AJ,] =3-p-p', (51b) 

i i 

Tr [7 • M±7o7 • fcA± 70] = Tr [7 • ^75 A±7 ' ^75 AJ,] = (1 + p • P')^^^fi^ (51c) 

It can be shown that the contribution from Eq. ()50ap to F is suppressed by an order g relative to that from Eq. (I50bp 
and wiU be ignored here. We neglect also the dependence of the gap function on the magnitude of the momentum, 
but keep the dependence on the energy and the momentum orientation. Then the integrals over p and p' can be 
carried out easily. Since we are only interested in the zero temperature, the Matsubara sum becomes an integral over 
the Euclidean energy. We find 

Tr[A.5n] = -%4 f I d'P I , ct^j^.m^^'^p') ^ [A(z.- 0) + A(^- 0)] (52) 

■^^'^ J J J J ^{i,^ + (f>^i,,p)){iy'2 + ^^i,',p')) 

where cos 6* = p ■ p' . Making use of the Nambu-Gorkov formalism in Eq. p3ll3^ . the rest terms of the condensate 
energy density Eq. ((32)) can be evaluated readily. 



^-Trln5-i - Trln^Q-i = J di, J d'p[\u\ - ^ly^ + (bHiy,p)\ (53) 

lTr(5-5 - 1) = (dJ dH^£t^ (54) 

The final expression of the condensation energy density reads 

f f f ,2^ f .i2Z,^r,.. ..' m ^{v.p)(j){v' .p') 



F = ~ ^-i^Jd:.Jd:.'Jd^pJ d^p'V{v~v',e) 
2^2 



^{v^^4>^{v,p)){v'^ + 4>^{v',p')) 



dv / d'p^==== (55) 



(27r)3 J J " ^ ^v-^ + ^-^{v,p) 
where V contains the contribution from both magnetic and electric giuons, i.e. 

V ^Di{v-v\e)+Dt{v-v',6) (56) 
The gap equation can be derived by minimizing F with respect to the gap function </)(:/, p), 

and we end up with 

0(.,p) fd.' [ d'p'V{.-.',9)— f^'''£\ (58) 

A consistent derivation of the gap equation up to the subleading order requires both the contribution from and 
that from the diagonal block of to be kept. The net result is to replace the first term inside the square root on 
RHS of (|58p by v'"^ /Z'^{v') with Z{v) the wave function renormalization of the normal phase. But it will not interfere 
with the angular dependence of the gap function to the subleading order as will be shown in the next section. 
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IV. THE INTEGRAL EQUATION FOR THE ANGULAR DEPENDENCE OF THE GAP 

Although the pairing strength are equal to the leading order of the QCD running coupling constant, similar to 
the case II of the toy model, the subleading terms fall off with an increasing J. This makes the solution to the gap 
equation (|58p highly nontrivial. In what follows, we shall isolate the energy dependence and the angle dependence of 
the pairing potential V. A differential equation with respect to the Matsubara energy will be derived from (ISS)) that 
fixes the gap function up to an arbitrary function of the angle. This function will be determined then by (|58p with 
a known function of the Matsubara energy. 

Proceeding with the partial wave analysis, we expand ¥{1^ — v' , 6) in scries of Legendre polynomials [l^]: 



V{v-v\e) = ^ln^^^(2J+ l)P;(cos0) + -L^(2J+ l)cjFj(cos(?) (59) 

^ ' J ^ J>0 

where ujc = "^"^'^t^'^ ^ and cj is given by Eq.([21). Using the completeness relation 

AT/ 9= 

^(2J+ l)Pj(cos0) = 47rj2(p_p') (60) 



J 

and the identity (proved in the Appendix A) 



j dH' p{2J + l)cjPj{p ■ p')f{p') = 2 j d^'l^^jl^ (61) 
with f(p) an arbitrary function of p, the gap equation (j58p becomes 



' In- ^ v-v , y 





1 



dx'- 



{v',x') <t>{v',x) 



\x — x'\ 



yjv'"^ + (p{v',x') ^Jv''^ + (^'^(y' ,x) 



where = 5^/(187r^), x = p - z with z a fixed spatial direction and a UV cutoff, ~ gjj, is introduced. In deriving 
([62|) . we have assumed that the gap depends on x only, so the integration over the azimuthal angle o f p' can be carried 
out explicitly. The gap equation (|62p can be further simplified by using the approximation of Son [17| 



In- ~ In- 63) 



with t/> = max(i^, J/'). It is convenient to introduce 



^ = \n—, a = ln^. (64) 



On writing (p — 0(^, x) and 



(65) 



the gap equation ([62|l becomes 



« 1 
^(c \ ccf^fc \ f ^^'^''^'^ , o f ^ , $(a,x')~$(a,x) 

m,x) = ^<S>{^,x) - J d^ ^ — +3 J dx i^^TT^Tj (66) 

a -1 

Taking the derivative of both sides with respect to ^, we find 

^-'^i^'^) (67) 
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which imphes the boundary condition 



- (68) 



as ^ ^ oo for aU x. Another derivative of (j67p yield the ordinary differential equation 

d^6 0^6 



= (69) 



which is universal for all x. It follows from Eq. (|66p that the gap equation is equivalent to 

1 

a^a,x)-(t>{a,x)+3 f ^^^4^^^^l^ = (70) 
J \x-x'\ 



The solution to ((69)) subject to the condition ((68)) contains an arbitrary function of x to be determined by ((70)) . No 
further approximation has been made up to now. 

The solution to the differential equation ()69p proceeds in the same way as that for a spherical gap. To the leading 
order, the equation can be approximated by a linear one, 

^ +5^0(6- 00^°) =0 (71) 

where b{x) is to be determined by the condition ii^^e'' = 1. Its solution that satisfies the boundary condition ((55)) 
and the continuity up to the first order derivative reads 

<j)o{x) cos g[b{x) - C], for ^ < 6(x), 
0(")(^,x) = <( (72) 
00 (x), for ^ > b{x). 



where 

h(T\ = In 

\M^)\ 



b{x) = Inrr^ (73) 



It follows from Eq. ((^ then that 



-^(j)o{x) sin g[b{x) - £,], for ^ < 6(x), 
H^,x) = { (74) 
0, for ^ > b{x). 



The angle dependent factor f{p) introduced in Eq.([2]) is defined by 



where Aq is the s-wave gap given by 



m - ^ = OH) (75) 



J-^ln^=0. (76) 



2 " Ao 



where the contribution from the wave- function renormalization is ignored. Up to the subleading order, the differential 
eauation ((M)) reads 

We find that 

0(1) (e, x) = 0(1) (e, x) + A{^, x)u{^, x) - B(e, x)v{t x), (78) 
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where u{(,,x) and v{^,x) are the two hnearly independent solutions to the Eq. ((7T 



and 



cos5[6(x) — for ^<b{x), 



1, 



for f > b{x). 



— sing[b{x) — for £_<b{x), 
for i > b{x). 



A{^,x)^g J de 0{b-e) 



v{S.\x)c^^°\£,\x), 



(79) 



(80) 



(81a) 



1 



B{^, x)^--g di' eib ~ a - 

At the point ^ = a, we have 



A{a,x) ^l + 0{g), B{a,x):^g\n2 

Therefore 

(j)^'^'^ (a, x) ~ 00 (2;) cos g{b — a) — .gln2 sin g{b — a) 

to the subleading order. Since g{b — a) = tt/2 + 0{g) according to Eq. ((76l) . we have 

<^{a,x) ^g(j)o{x) + 0{g), 

and 

0(1) (a, x) ~ 00 (x) [| - g{b ~a)~ .gln2] + 0(5) 
Substituting Eas. (|84)) and ((SS]) into Ea. (|70|) . we obtain the gap equation to the subleading order 



TT _ _ UJc 

2 ^''|0o(^)| 



5ln2 



0o(x) + / dx ^— = 



Then the integral equation for f{x) 



fix)ln\f{x)\-3 I dx'l^^^^LJ^^^o 

\x — x'\ 



(81b) 

(82) 
(83) 
(84) 
(85) 

(86) 



(87) 



follows from ([76|) . 

Few comments are in order: 1) The spherical gap, f{x) = 1 is a trivial solution to Eq. (|87p and there is no angular 
momentum mixing. 2) The "spin-1" gap, carried over from the polar phase of He'^, f{x) oc x, fails to satisfy this 
equation. 3) Ea. (j87p conserves the parity. In another word, its solution can be either an even or an odd function 
of X. 4) If the wavefunction renormalization is restored, there will be an additional subleading term on RHS of (j77[) 
and an additional subleading term on RHS of B of eq. ([82| . This term, when substitute into Eq. ([70| . will cancel the 
corresponding contribution to Ao 



Ao = ^e-^Ti°) 
leaving the integral equation ([57)1 intact. 



20487rV _3^_£+± 
— — e ^ 9 

Njg^ 



(88) 
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V. THE NUMERICAL RESULTS OF THE ANGULAR DEPENDENCE 



The solution to the integral equation Eq. (j87[) can be obtained from a variational principle. Upon substitution 
of Eg. ([7^ with (t)o{x) = Ao/(a;) into Eg. ([55)1 . the condensate energy density becomes a functional of /(details in 
Appendix B), i. e. 



where 



T[f{x)] = I dxf\x) ln|/(a:)|-- 



2 J J \x-x'\ 

-1 -1 



(89) 



2 / dxrix) \n\fix)\- 



1 



3 J dx J dx I 





,([f{x)-f{x')Y , [f{x) + f{x')\ 



\x — x'\ 



X -\- x' 



(90) 



with the last equality following from the odd parity of f{x), i. e. f{^x) — ~-f{x). Readers may easily verify that the 
variational minimum of Eg. ((90)) does solve Eg. (|87)) . 

Before the numerical solution, we consider a trial function 



f{x) = cx 

and substitute it into the target functional ()90p. The minimization yields 

c = e"^ = e-^+^ ~ 3.459 x lO^^^ 

at which 

T[f] ~ -3.989 X IQ-^ 



(91) 
(92) 
(93) 



The trial function ((9T)) is what people carried over from the polar phase of He^. The "-6" of the exponent of ((92)) 
comes from the pairing strength of the p-wave and the " i" stems from the logarithm of ((57)) . The latter contribution 
was reported in 13]. The trial function ()9ip with p-wave alone is not optimal. The free energy will be lowered further 
by including higher partial waves of odd J as we shall see. 

To find the variational minimum, we discretize the integral of Eg. ((90|) by dividing the domain x £ (0,1) into 
N{>> 1) intervals with 



Xj^{j + l)Ax, i = 0,l,2,---,iV-l 



where Ax = 1/N. We have then !F — limAr^oo with 
.fAr-2Ax^/2(ln/ 



u ..... L \Xi Xj Xi -\- Xj 



(94) 



(95) 



where we have dropped the limit x' ^ x of the first term inside the curly bracket of Eq. ((90)) . is a function of N 
variables. The stationary condition 







yields 



fjlnf,+3Ax\— + J2^ ^ 



1 



Xi Xj \ Xi ~\- X j 



1 



\xi Xj I Xi -\- Xj 



/. = 



(96) 



(97) 



which is a discrete version of Eq. ((57)1 . Regarding /^'s as given, the equation for fj is of the form 

(In/, + a)/, -6 = 



13 





T 1 

J — 1 


T q 

J — O 


J — 


bj 


3.413 X 10"^ 


-2.328 X 10"^ 


7.409 X lO"'^ 



TABLE I: The first three expansion coefficients of the gap function according to Legendre polynomials. 
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FIG. 2: The angular dependence of the gap function with angular momentum mixing. The dashed line and the solid line are 
the initial configuration and the final numerical results respectively. 

with a and b positive. It has one and only one solution for fj > 0. 

We start with the trial function (PT|) . fj = e~^Xj as an initial configuration and update each fj by solving Eq. ([98|) . 
This way we lower the value of the target functional in each step and approach the solution to ([87|) eventually. The 
process converges rapidly and our numerical solution to (j87p is shown as solid line in Fig[21 which depart from the 
trial function (dashed line) slightly. We find the minimum value of the target functional 

T[f] ~ -4.130 X 10"^. (99) 

which drops from ([93|) by 3.5 percent. 

It is instructive to examine the angular momentum contents of our solution in the partial wave expansion 

.f{x)^ bjPj{x). (100) 

J=odd 

The coefficients of the first three partial waves, J = 1, 3, 5, calculated by substituting the numerical solution into the 
formula 

1 

2 7+1 f 

bj = ^^ J dxf{x)Pj{x) (101) 
-1 

are displayed in Table I. While the gap function contains all partial waves of odd J, the component of J = 1 is the 
biggest. This is anticipated because the pairing strength of the all partial waves are equal in leading order but fall off 
with an increasing J in the subleading order as is shown in the partial wave expansion ([55)) . 
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VI. CONCLUDING REMARKS 



In summary we have explored the angular dependence of the gap function for a non-spherical pairing of CSC. 
Because of the equal strength of the pairing potential mediated by one-gluon exchange for all partial waves to the 
leading order of QCD running coupling constant and the nonlinearity of the gap equation, a non-spherical gap function 
cannot be restricted to one angular momentum channel only. Other multipoles are bound to show up, which renders 
the angular dependence of the gap nontrivial. On the other hand, the pairing strength to the subleading order 
decreases with increasing angular momentum J . The mixing effect will not be as big as that in the soluble toy model 
we introduced for the purpose of clarification. 

For the single flavor CSC, we worked out the angular momentum mixing effect explicitly for the gap function 
with zero azimuthal quantum number at zero temperature. An nonlinear integral equation for the nontrivial angular 
dependence was derived and its solution was obtained numerically. The gap function in this case reads 



where Aq is given by Eq. ([88|) and /(p • z) is plotted in FiglS] 

The drop of the free energy of the modified polar phase by the mixing, however, is numerically small. The magnitude 
of its condensation energy is smaller than that of the CSL phase by a factor of 1.48 instead of the factor 1.54 reported 
in [l3]. The CSL phase remains stable. In this sense our results at the moment is of theoretical values only. There 
are many other candidate pairing states between quarks of the same flavor [l3, [l^ . Among them are the states with 
a nonzero azimuthal quantum number and the pairing between quarks of opposite helicities. The former is analogous 
to the A phase of He'^ and may be present in a compact star with a strong magnetic field. The pairing force in the 
unequal-helicity channel is stronger [l^ . fisl . [2^ . The angular momentum mixing effect is generic in all nonspherical 
pairing states and the integral equation ([57)) can be readily generalized to these cases. There may be phenomenological 
implications of the angular momentum mixing. A systematic survey of the angular momentum mixing effect in all 
"spin-1" CSC states covered in 13] will be reported in another paper. 

Another place where the angular momentum mixing shows up is the CSC-LOFF state in the presence of Fermi 
momentum mismatch. It has been speculated [28| that the forward singularity will increase the upper limit of the 
mismatch value that supports a LOFF pairing. The new threshold was found in ^29j , motivated by the nearly equal 
pairing strength of all partial wave channels. The same mechanism works for the gap equation of LOFF pairing. Its 
free energy will be lowered by the angular momentum mixing and the lower edge of the LOFF window is expected to 
be shifted to a lower value of the mismatch parameter. 
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APPENDIX A: THE DERIVATION OF EQUATION (f6T|) 



The integral formula of cj is [2! 



1 




(Al) 



-1 



It is convenient to introduce 



1 



1 




(A2) 
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where e(> 0) is an infinitesimal quantity. We have \im^_^Q+ Cj = cj. For the first term on RHS, we expand 



according to 



2J+1 f Pj{x 
aj = — : — / ax 



Therefore Eq. (|A2|) reads 



2 ./ l~x + e 
1 



e 2 2 + e 

Cj — ^ aj — m- 



2J+1 2 
Evaluating the summation in Eq. (j6ip is straightforward now 



J2(.2l + ^yjPj{p-p') = ^2a,7Pj(p-p')-^(2J+l)Pj(j5.p')ln 



J=i 



J=i 



1 — p • p' + e 



— 47r5^(p — p')ln 



2 + e 



Then for an arbitrary function /(p), 



d'p'f2{2J+l)c:jPj{p-p')f{p') = 2 / d^T '^^^'^ 
J=l 



47r/(j5)ln 



2 + e 



2 / 



J{p')-f{p) 
I ~ p ■ p' + e 



(A3) 



(A4) 



(A5) 



(A6) 



1 — p ■ p' + e 



47r/(p)ln 



2 + e 



2 / 



2V/(P')-/(P) 



1 — p • + e 



(A7) 



The Eq. ((6T|) is obtained by taking the limit e ^ 0^ 



APPENDIX B: THE CONDENSATION ENERGY DENSITY WITH THE ANGULAR MOMENTUM 

MIXING 

In this appendix, we shall derive the expression Eq. (|55p ) of the condensation energy density with the angular 
momentum mixing. Substituting Eq. (|59p into the first term of Eq. (|55p , we find 



327r4 
1 



d.'J d. J dH' j ^'p[6^i"^;^^(2Z + i)PK75-p') 



2/i2 

167r3 



^(2/ + 1)qPKp-p') 



4){v,p)4>{v' ,p') 



1=1 



dv' j dv j d^pln 



<^c (j){v,p)(j){v' ,p) 



-jd.'jd.jdH'jSp^ 



I 



(t){v,p) 



Because of the eveness oi (j){v,p) in v, we have 

Wo Wo 



f Id^pl d.' I In'^' 



4-^-3 



(j){v,p)(j){v' ,p) 



{v^u')]] (Bl) 





Wo Wo 



(l){v,p) (j>{iy,p' 



(v ^ v') I (B2) 
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where the approximation (|63p has been appHed to the forward logarithm. For the gap function of zero azimuthal 
quantum number, p) depends only on x = p ■ z. We find that 



1 C30 oo 



1 1 



- 1 a a 



dx 



-1 -1 



[$(a,a;') - $(a,a;)]^ 



(B3) 



where, $(^,a;) has been defined in Ea. (p5)) and ^ and a have been defined in l|64p . The integral over ^' followed by 
the integral by part over ^ leads to 



1 C30 



1 1 



Fi = ^-^^{ - a J dx<^^{a,x) - J dx J d^^^{^,x) + ^ J dx J dx 



,[a>2(a,a;')-$2(a,a:) 



■} (B4) 



Making use of Eq. ([65| and ([85|) , we have 



d^^'iC, x) = -r^oix) --g{b-a)- g\n2 



f / dS, 



c^\^.x) 



I _|_ 0^({,a:) ^2f 



(B5) 



and thus 



a^\a, x)- j d^^\^, x) = g4>l{x) {\ ~ gb ~ glni) ~ f j dy 



4>^{v,x) 



+ (/)2(z/, x) 



(B6) 



Substituting (/>o(x) = Ao/(a;) into (|B4|) . we obtain that 



27r2 



1 1 



dxf{x)\n\f{x)\ 



dx i dx 



,[f{x)-f{x')] 



2^2 



-1 -1 



^i^2 ^ 02 (-j^^ 



(B7) 



Then the condensate energy density with the angular momentum mixing reads 



F = Fi 



47r3 



dv j d p 



27r2 



dT/2(x)[ln|/(a:)|-l 



1 1 



— I dx I dx 



, [fi x) - fix') 
\x — x'\ 



(B8) 



The minimization of this free energy give rise to Ea.([57|) of the text. 
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